Multiple plasmon resonances in naturally-occurring multiwall nanotubes: infrared 

spectra of chrysotile asbestos 
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Chrysotile asbestos is formed by densely packed bundles of multiwall hollow nanotubes. Each 
wall in the nanotubes is a cylindrically wrapped layer of Mg3Si205(OH)4. We show by experiment 
and theory that the infrared spectrum of chrysotile presents multiple plasmon resonances in the 
Si-0 stretching bands. These collective charge excitations are universal features of the nanotubes 
that are obtained by cylindrically wrapping an anisotropic material. The multiple plasmons can be 
observed if the width of the resonances is sufficiently small as in chrysotile. 

PACS numbers: 78.67.Ch, 63.22.+m, 73.20.Mf, 91.60.-x 



The interaction of nanostructures with electromagnetic 
waves has recently received an increasing attention, both 
because of fundamental and technological issues. The 
frequency dependence of this interaction is governed by 
different phenomena according to the characteristic size 
of the nanostructure template. The regime dominated by 
interference and diffraction occurs for particle sizes close 
to the wavelength of light as, e.g., in photonic crystals or 
in opals. The regime dominated by quantum mechanics 
occurs for sizes close to or smaller than a few nanome- 
ters such as in quantum dots or nanocrystals. Between 
these two regimes, a third regime occurs where the fre- 
quency dependence is dictated by collective charge exci- 
tations, i.e. by confined plasmons p], 0, 0, [| [|. In 
this regime, an interesting phenomenon occurs in multi- 
wall nanotubes that are obtained by wrapping around 
cylinders layers of an anisotropic material, as in car- 
bon jl], [|, H or WS 2 H nanotubes. In particular, the 
macroscopic dielectric tensor, *e(u>), of the nanotubes can 
present multiple confined plasmon resonances, that corre- 
spond to a single resonance in the dielectric tensor of the 
same material with flat planar layers. Such phenomenon 
has been observed in theoretical calculations ||. Exper- 
imental measurements exist in the energy range of elec- 
tronic excitations (ui ~ 5-25 eV) [|, [|, ||, ||. However, in 
these experiments, the multipeak spectrum exhibits only 
one broad resonance because of the large width of the 
electronic resonances of the corresponding planar mate- 
rials m. 

By contrast, the resonances oi*e(uj) related to the ionic 
vibrations, in the infrared (IR) range (lu ~ 0.1 eV), are 
ideal candidates for the detection of multipeak spectra 
produced by collective charge excitations. In fact, the 
ionic resonances are much sharper than those related to 
the electronic excitations. Thus, in ionic materials with 
intense IR activity, multiple plasmons could be expected. 
Here, we show that the IR spectra of chrysotile asbestos 
present, indeed, multiple plasmon peaks related to its 



FIG. 1: Multiwall nanotubes of chrysotile, arranged as a 
hexagonal close packed array. The cylinders are infinitely 
long in the z direction. Inset: Structure of lizardite showing 
the stacking of two layers. 



peculiar tubular nanostucture. 

Chrysotile asbestos is an example of natural mul- 
tiwall hollow nanotubes with Mg3Si20s(OH)4 compo- 
sition, belonging to the group of serpentine minerals. 
These minerals are formed by stacking layers, each con- 
taining a pseudo-hexagonal silica sheet of corner-shared 
Si04 units linked to a trioctahedral sheet of edge sharing 
Mg0 2 (OH) 4 octahedra (Fig. 0). The OH groups located 
at the top of the trioctahedral sheet are H-bonded with 
the O of the basal plane of the next layer. These minerals 
present cylindrically wrapped layers in chrysotile, corru- 
gated layers in antigorite, and flat layers in lizardite @, |). 
The structural relation between chrysotile and lizardite 
is the same as that between C nanotubes and graphite. 
Chrysotile nanotubes present outer and inner diameters 
of 240 - 600 A and 50 - 100 A, respectively (Fig. |). The 
tubes are arranged as densely packed bundles of parallel 
fibers, forming veins up to several centimeters thick. 

We measure the IR spectra, in transmission and at- 
tenuated total reflectance (ATR), of a chrysotile sam- 
ple from the asbestos mines of Salt River Canyon, Ari- 



FIG. 2: Transmission electron micrography of the investi- 
gated chrysotile sample, perpendicular to the tube axes. In- 
set: selected area electron diffraction pattern showing the par- 
allel texture of chrysotile tubes seen along [100]. 




FIG. 3: Transmission powder IR absorption spectra of 
chrysotile and lizardite in KBr pellets (absorbance units). 
Spectra are recorded at room temperature using an FT-IR 
spectrometer Nicolet Magna 560 with a resolution of 2 cm . 
z and p denote polarizations parallel and perpendicular to the 
tube axis, respectively. 



zona. The sample, characterized by transmission elec- 
tron microscopy and electron microprobe, consist of pure 
chrysotile tubes in parallel texture (Fig. |J). The IR 
powder spectrum of chrysotile, compared to that of the 
flat-layered variety lizardite, displays a strong absorption 
band at 1018 cm -1 in the range of Si-0 stretching vibra- 
tions (Fig. |) (see also |, ^ [Q]> (L2|). The shoulder 
observed at 1018 cm -1 in the lizardite spectrum is in- 
deed related to chrysotile and polygonal serpentine con- 
taminations present in low amount in natural lizardite 
samples [[ll, 13 . In the IR powder spectrum of lizardite, 
the band at 1084 cm -1 (Li, Tab. |) is related to the 
out-of-plane (perpendicular to the layers) Si-0 stretching 



FIG. 4: Attenuated total reflectance IR spectra of chrysotile. 
Spectra in absorbance units are recorded with an ATR Thun- 
derdome device equipped with an AgBr polarizer. The orien- 
tation of the chrysotile aggregate with respect to the incident 
IR wave is indicated; Chr: chrysotile; Ge: germanium crystal. 



mode, whereas that at 951 cm" 1 (L 2 , Tab. |) is related to 
the two degenerate in-plane (parallel to the layers) Si-0 
stretching modes |L3|. The bands at 961 and 1018 cm" 1 
of chrysotile were previously related to the lift of this de- 
generacy by the bending of the tetrahedral sheets |To[ . 
We better constrain the dependence of the IR spectrum 
as a function of the polarization of the incident wave, by 
recording the ATR IR spectra of an oriented chrysotile 
aggregate. In this experiment, a thick oriented sample of 
chrysotile is pressed onto a Ge crystal. The spectrum of 
the IR light propagating in the Ge crystal and reflected at 
the chrysotile/Ge interface is recorded. The IR light has 
an incidence angle of 45°. The propagation and electric 
field vectors of the incident wave are set either parallel or 
perpendicular to the fiber axis and to the incidence plane, 
respectively (Fig. ||). The ATR spectrum recorded with 
the electric field vector parallel to the fiber axis displays 
only one band, whereas the three other spectra display 
three bands. 

Since the size of the nanostructure template is much 
smaller than the wavelength of light, the IR spectrum 
is fully defined by the macroscopic dielectric tensor of 
chrysotile < i*c(w). To model the spectra, we consider, 
for simplicity, a nano-structured material composed of a 
close- packed array of mono-dispersed nanotubes (Fig. |l|). 
An accurate effective medium theory has been derived for 
such systems pL However, here we compute *ec( w ) from 
the polarisability of a single tube using a simple Clausius- 
Mossotti model, given the approximations made on the 
size distribution and on the arrangement of the tubes. 

The polarisability a(uj) is calculated from the macro- 
scopic dielectric tensor of the planar material, lizardite. 
The dielectric properties of an isolated tube in vacuum 
can be derived by assuming that the chrysotile tube is a 
dielectric continuum, locally identical to lizardite. Thus 
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the layer curvature is accounted for only at the macro- 
scopic scale, as in Refs. [0, |], ||. The low-frequency 
dielectric tensor of lizardite, can be written in 

Cartesian coordinates as: 

7 l {uj) = € ± (lu)u x u x + ej_(uj)u y u y + £||(£J)U Z U Z , (1) 

where u z corresponds to the c axis of the lizardite hexag- 
onal structure, and u> is the frequency of the IR wave. In 
Ref. (l3), we have obtained the theoretical < ?l(w) from 
the dynamical matrix and the effective charge tensors. 
These quantities were calculated from first-principles us- 
ing density functional perturbation theory p| . In the 
present study, we use £i(w) of Ref. [jl3| calculated with a 
damping coefficient of 16 cm" 1 . However, to better fit the 
experimental IR spectrum of lizardite, we shift the trans- 
verse optical frequency of the two degenerate in-plane 
Si-0 stretching modes from the theoretical value of 915 
cm" 1 |lj| to the experimental value of 951 cm" 1 . The 
resulting theoretical IR powder spectrum of lizardite, ob- 
tained following M, is reported in Fig. ft. 

The microscopic dielectric tensor of a chrysotile tube 



ec(u), is then obtained using an effective medium ap- 
proach of Clausius-Mossotti type j| : 



'e, 1 c{<P,r, w) = e±(uj)u z u z + e||(w)u,.u r + e ± (w)u v u ¥ 



(2) 



for Ri < r < i?2, and e^cif, t, u>) = 1 otherwise. Here 
u z , u r and u^, are now the unitary basis vectors of cylin- 
drical coordinates and R\ and R2 are the inner and outer 
radii of the tube. 

The polarizability of a single tube is obtained from the 
the electric field induced by a quasi-static external field, 
E oxt . The component of the electric field parallel to the 
cylindrical axis is homogeneous. The perpendicular com- 
ponent of the electric field is E(y, r, u>) — —WV((p,r,u)). 
The electric field corresponding to a perpendicular exter- 
nal field E oxt is derived from the first Maxwell equation: 



V • \?»G(.<p,r,w)E(<p,r,u)] = 0, 



(3) 



with the boundary condition E(y, 00, to) — E cxt . The 
electrostatic potential solution of Eq. (|^) has the form: 

V(ip, r, u>) = air cos p, (4) 
V{ip,r,u) = [a 2 r A{u}) + 6 2 r~ AM ] cos <p, (5) 
V(p, r,u>) = (- E cxt r + b 3 r~ 1 ) cos tp, (6) 

for r < Ri, Ri < r < R 2 and r > R 2 , respec- 
tively. Here A(w) = y/e±(u}) /e\\(u)) and the scalar 
quantities 01,02,62 and 63 are determined by the con- 
tinuity of V(ip, r, to) and of the radial component of 
eVcfVi r i w )E(i J 9, r, uj) at the inner and outer surfaces of 
the tube. From these expressions, the transverse polar- 
izability of the single tube a p (oj) can be derived (Eq. (4) 
of pi). The macroscopic dielectric tensor of chrysotile, 



ec,*M - f'e±(u) + (1 - /'), 



l-/op(u;)/2' 



(7) 
(8) 



where ec,z(^>) and ec, p {oS) correspond to the components 
of ^cfw) parallel and perpendicular to the cylindrical 
axis, respectively, / is the volume fraction occupied by 
the tubes, / ~ 0.9 for a close packed arrangement of the 
tubes and /' = /[l — {Ri/R 2 ) 2 ]- From these equations, 
it can be shown that a single resonance in the dielec- 
tric tensor of the flat planar material, *ei,{u)), can pro- 
duce multiple resonances in a p (w), and therefore in the 
macroscopic dielectric constant ec&ioj). This behavior 
has been observed in Fig. 6 of Ref. ||. The frequency for 
which confined plasmons occur, can be obtained as the 
non-zero solutions of Eqs. (H), (||) and (||) with E cxt = 
0. In the limit of a vanishing damping coefficient, these 
resonances occur when D(u>) = — zA(w) is real and satis- 
fies: 

D(w)log(R 2 /Ri) = -2atan[e||(w)£)(cj)] +kn, (9) 

where k is an integer. This equation admits an infinity 
of solutions because of the divergence of ^(w) at a res- 
onance with zero damping. In practice, the finite width 
of the resonances of *e l(uj) leads to a finite number of 
resonances of *ec(o;), increasing with a decreasing ratio 
R\/R 2 . 

To calculate the ATR spectra, we compute the reflec- 
tion coefficient of an electromagnetic wave incident on 
a flat interface between the chrysotile medium, with di- 
electric tensor ec{oS), and a Ge medium, with isotropic 
dielectric tensor equal to 16. The theoretical spectrum 
in absorbance unit is obtained as minus the logarithm of 
the reflection coefficient (Fig. ^). The major features and 
the changes of the spectra as a function of the polariza- 
tion of the incident wave are very well reproduced by our 
calculations, with a ratio R\jR 2 = 0.1, comparable with 
that observed by TEM (Fig. ||). The spectrum obtained 
with the electric field parallel to the fiber axis (Fig. |J, 
1 st panel) is described by ec, z (w) and displays a single 
band (C4, Tab. |). The spectra obtained with the electric 
field perpendicular to the fiber axis (Fig. |J, 2 nd and 3 rd 
panels) are described by ec,p(kO and display three reso- 
nances (C3, C 2 , and Ci, Tab. |). Finally, the spectrum 
obtained with the geometry of Fig. |], last panel, contains 
resonances related to both ec, p (u>) and cc,z(u). 

The theoretical powder IR transmission spectrum of 
chrysotile is computed using as in Ref. jljj]. In 

this case, we consider that micrometric chrysotile par- 
ticles are elongated along the z axis and inserted in a 
homogeneous KBr medium. An excellent agreement be- 
tween the theoretical and the experimental spectrum is 
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FIG. 5: Dissipated power density in the nanotube at the 
resonances C2 and C3. The darker gray indicates larger dis- 
sipation. The arrow shows the direction of E cxt . 

obtained (Fig. ||). In particular, the calculation repro- 
duces the presence and the position of the strong absorp- 
tion band at 1018 cm -1 , characteristic of chrysotile. By 
an analysis of the theoretical spectra, we can correlate 
the peaks observed in ATR to those observed in trans- 
mission. In particular, we find that the peaks at 1084 
and 1018 cm -1 in transmission are associated to the res- 
onances Ci and C2, whereas the resonance at 961 cm -1 
is mainly associated to C4, and in a minor proportion 
to C3. The frequency shifts and intensity changes of the 
bands between the ATR and the transmission spectra 
arise from the different geometry of the two experimen- 
tal setups, and are reproduced by our calculations. 

Our theoretical investigation of the IR spectrum of 
chrysotile also provides an unambiguous assignment of 
the resonances with respect to the vibrational modes of 
the layers, Tab. |. The out-of-plane Si-0 stretching mode 
leads to Ci. The in-plane Si-0 stretching mode polarized 
parallel to the cylindrical axis of the tube leads to C4. 
The other in-plane Si-0 stretching mode, perpendicular 
to the cylindrical axis, leads to the two resonances C2 
and C3, that correspond to the solutions of Eq. (||) for 
k—1 and k=2, respectively. Therefore, the resonances 
C2 and C3 are related to two distinct localized plasmons 
which originate from the same resonance in the dielec- 
tric tensor of lizardite. We study the spatial extension 
of these two plasmons by looking at the spatial distri- 
bution of the dissipated power density, W(<p,r,uj) [ El| . 
In Fig. |^, we report W(cp,r,u>) for an isolated tube for 
the resonances C2 and C3. The resonance C2 is related 
to a surface plasmon localized at the internal surface of 
the nanotube. The resonance C3 is distributed over the 
whole tube, and presents two nodes close to the inner 
surface (Fig. ||). Note that, in our calculations, the two 
in-plane Si-0 stretching modes of the layered local crys- 
tal structure are still degenerate. Therefore, the splitting 
of the Si-0 stretching absorptions bands in the IR spec- 
tra of chrysotile is not related to a change in the micro- 
scopic structure of the layer induced by the bending of 
the tetrahedral sheets, as previously proposed flQfl . 

In conclusion we have shown, by experiment and the- 
ory, that multiple plasmons are observed in the IR spec- 
tra of chrysotile nanotubes. In particular a single in plane 
Si-0 stretching mode of the Mg 3 Si205(OH)4 layer pro- 
duces two resonances in the polarisability of chrysotile. 



TABLE I: Theoretical transverse optical (TO) modes and 
corresponding resonances of the IR spectra. Experimental 
values in parenthesis. With || and _L we indicate the Si-O 
modes parallel and perpendicular to the Mg3Si205(OH)4 lay- 
ers, and with u r , u v , and u z , the dominant cylindrical polar- 
ization. 





polarization TO 


label transmission ATR 


Lizardite 




1045 


Li 


1105 (1084) - 




_L 


951 


L 2 


951 (951) 


Chrysotile 


II Ur- 


1045 


Ci 


1077 (1084) 1071 (1075) 






951 


c 2 


1022 (1018) 992 (995) 






951 


c 3 


972 955 (955) 




_L u z 


951 


c 4 


951 (961) 936 (938) 



We have predicted that the number of resonances in- 
creases with a decreasing width of the vibrational modes 
and/or a decreasing ratio between the inner and outer 
radii of the nanotubes. Thus we anticipate that more 
resonances should be detected in chrysotile samples or 
other nanotube systems which optimize these parame- 
ters. Finally the IR plasmons resonances could become a 
tool to monitor the geometrical parameters of nanotube 
based materials. 

Calculations were performed at the IDRIS, and TEM 
work at the French INSU National Facility, CRMC2- 
CNRS, Marseille. This is IPGP contribution #0000. 
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